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We present a general algorithm which permits to construct solutions in string cosmology for 
heterotic and type-IIB superstrings in four dimensions. Using a chain of transformations 
applied in sequence: conformal, T-duality and SL(2, R) rotations, along with the usual 
generating techniques associated to Geroch transformations in Einstein frame, we obtain 
solutions with all relevant low-energy remnants of the string theory To exemplify our 
algorithm we present an inhomogeneous string cosmology with S s topology of spatial 
sections, discuss some properties of the solution and point out some subtleties involved in 
the concept of homogeneity and isotropy in string cosmology. 
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1. Introduction 



Cosmological implications of string theory are receiving quite a broad attention these 
days (for a necessarily incomplete list of recent and not- so-recent references see @ @ 0)- 
This is not by chance since, given the lack of traces of stringy effects in "low energy" 
particle physics, the cosmological scenario appears as a perfect arena for the search of 
string blueprints in the early Universe. 

In spite of the impressive progress in the study of string theory beyond the perturbative 
regime (for a review see, for example, ||), we still lack a full-fledged non-perturbative 
formulation which allows a description of the early Universe at Planck time, where the 
classical concepts of space and time cannot be used any longer. In the meantime, one 
is bound to study classical cosmology using the low-energy effective action induced by 
string theory, which generalizes General Relativity by including other massless fields. As 
a consequence, in string cosmology one continues to ask essentially the same questions 
earlier formulated in the framework of Einstein Relativity, but now posed in the presence 
of those extra degrees of freedom. The central puzzle of theoretical cosmology remains as 
well the same, and is concerned mainly with the question of why the present Universe looks 
as isotropic and homogeneous at large scales. The desire to answer this last question leads 
one to consider the initial conditions and the evolution of the Universe at the beginning of 
expansiorJI. 

Apart from a few exceptions most of the work in string cosmology is focused 

on homogeneous solutions. However, since in General Relativity the generic solution near 
the cosmological singularity is neither isotropic nor homogeneous, one is led to study 
cosmological models with less symmetry, leaving aside not only the principle of isotropy but 
homogeneity as well. The ideas to study such a general behavior were inspired by Landau 
|J and worked out in details in classical papers by Belinski, Lifshitz and Khalatnikovil 
(BLK) [p!D| . The techniques due to BLK are strongly rooted in physical intuition, but lack 
however rigor when approximations are used in a highly non-linear regime. Although the 

We will not dwell here on different approaches and views as to how to solve this problem, but 
just mention some central ideas such as chaotic cosmology program H, C-gravitational entropy 
@ and the inflationary scenario || . 

5 In string cosmology the BLK approach has been recently discussed in the context of the 
Pre-Big Bang scenario by Veneziano [|J]. 
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conclusions of this analysis are probably correct at large, one must necessarily confront 
them against the analytic behavior of exact solutions. 

Because of the mathematical complexity of the generic inhomogeneous models, one 
usually deals with solutions in which homogeneity is broken in only one direction. The 
space-times obtained this way are usually referred to as G2 or Einstein-Rosen cosmologies 
TUl ||12|1 1|13| and are thought to provide the leading approximation to a general solution near 



the initial singularity [|14j]. It is important to mention here the relation of the Einstein- 



Rosen cosmologies to the interesting problem of scattering of plane waves in General Rel- 
ativity [ D)| [|16| . This relation, considered first in [|I7] , consists just in inverting the arrow 



of time of the vacuum, or vacuum with massless fields, inhomogeneous solutions. After 
all, the approach to the focusing two-surface in the collision of self-gravitating plane waves 
(not necessarily gravitational) is just the time inverse of the behavior near the singularity 
in a G2 cosmology. This being so, the study of the behavior of string induced G2 cosmolo- 
gies may lead to a better understanding of such a highly non-linear process as plane wave 
scattering. 

The main purpose of this paper is to describe a sufficiently general algorithm which 
permits the construction of new solutions of string cosmology starting from vacuum so- 
lutions of the Einstein field equations. Vacuum G 2 space-times will be used as the basic 
building blocks of our string solutions. Before going any further in the description of our 
generating technique, it will be useful to briefly mention some facts about the low-energy 
effective theory of the different superstring models and their symmetries. 

1.1. String low-energy effective actions 

On general grounds, the massless bosonic sector of superstrings includes the gravita- 
tional field G^v, the dilaton (p, with vacuum expectation value determining the string cou- 
pling constant, and the antisymmetric rank-two tensor B$ , besides other fields depending 
on the particular superstring model under study. The lowest order effective action for the 
massless fields can be written as 



'eff = — Aor / d D xVGe- 2 * R + 4(dcP) 2 - ±-{H^) 2 
(a ) 2 J I J-^ 



+ Smd, (LI) 



where H^ 1 ' = dB^ is the field strength associated with the NS-NS two-form and S^d is 
a model-dependent part which includes other massless degrees of freedom. When D < 10 
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some of these massless fields correspond to gauge and moduli fields associated with the 
specific compact ificat ion chosen, and the dilaton <f> appearing in is related to the 

ten-dimensional dilaton </>io by 2<p = 2<pio — logFio-D, where Viq-d is the volume of the 
internal manifold measured in units of \fa' ' . 

We restrict our attention to generic degrees of freedom, leaving aside the internal 
components of the ten-dimensional fields. In the heterotic string case, S^d contains the 
Yang-Mills action for the background gauge fields A®, which we set to zero in the following. 
In the case of the model-dependent part of the type-IIB superstring things are slightly more 
involved; among the massless degrees of freedom in the R-R sector we find, along with a 

(2) 

pseudo-scalar x (the axion) and a rank-two antisymmetric tensor B^J , a rank four self- 
dual form A sd urjX . The presence of this self-dual form spoils the covariance of the effective 
action for the massless fields, since there is no way of imposing the self-duality condition 
in a generally covariant way. But, if we set A sd to zero, we can write a covariant action 
for the remaining fields with 

-fUB 



QlLti 



(a'Y 

Notice that the R-R fields do not couple directly to the dilaton. Thus, the lower dimensional 

(2) 

(D < 10) R-R fields x an d B^J are obtained from the ten-dimensional ones through 
X = V^DXio and B& = ^/V^b[ 2} . 

By combining the dilaton and the axion of the ten-dimensional type-IIB superstring 
into a single complex field A = xiQ + ie~^ 10 , it is possible to check that the bosonic 
effective action is invariant under the £X(2,R) transformationi A — > (aX + b)/(c\ + cf), 
G^v — > |cA + d\G^ v ||18|1 . Writing this transformation in terms of four-dimensional fields 
(D = 4) we find that the action is invariant under the following field redefinitions (cf. J11 

/ bd + ac e~ 2 ^ 4 ac%| + (ad + bc)xA 

XA 77 . tTt . T OA, 1 1 



[(cxa + d) 2 + c 2 e~ 2 ^] i [(cxa + d) 2 + c 2 e~ 2 ^} \ 



[(cxa + d) 2 + c 2 e" 2 ^]i 
HW =dHM-cHW 

H W = [(c X A + d) 2 + c 2 e- 2 ^}H-bH^+aH^) 
G'^ = [(c X A + d) 2 + c 2 e- 2 ^G^ 



C ( 1 ") f 2") 

The two antisymmetric tensors B^J and B^J also transform as a doublet under SL(2, R). 
Notice that here we are using the string frame in which the metric is not invariant. 
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Comparing with the corresponding transformations for the ten-dimensional fields we find 
that the four- dimensional ones acquire an "anomalous weight" due to the fact that Vq, as 
measured in the string frame, does transform under SL(2, R). 

The algorithm to be described in the following uses as input diagonal vacuum solutions 
to Einstein equations with two commuting space-like Killing vectors. Infinite dimensional 
families of solutions of this kind are known [11] [ II| [17| [|13 . After transforming these diag- 
onal vacuum solutions to generate off-diagonal terms in the metric, we generalize them to 
include a minimally coupled massless scalar field. A conformal rescaling of the Einstein 
metric gives solutions to four-dimensional dilaton gravity and a T-duality transformation 
generates an antisymmetric rank-two tensor Bj^J from the off-diagonal components of the 
metric. This leaves us with string cosmology solutions of both the heterotic and type-IIB 
superstring. In the latter case, we still can generate non-trivial background R-R fields 
by performing a SL{2, R) rotation of the solution (cf. [|J). In the next section, we will 
detail the algorithm and in Section 3 a concrete application will be worked out. Finally, 
in Section 4 we will summarize our conclusions. 



2. String cosmology from General Relativity 

After having outlined the key steps to be followed in the construction of solutions, 
we now turn to the generating algorithm itself. In order not to get too involved with the 
details we shall ommit some technicalities, referring in cases to the original papers. We 
present here only those details and expressions necessary to get the final solution and to 
keep the presentation self-consistent. As mentioned above our starting point is a globally 
diagonalizable G<i line element which may be put into the following convenient form 

ds 2 = e f(t > z) (-dt 2 + dz 2 )+-f ab (t 1 z)dx a dx\ a, b =1,2. (2.1) 

The local behavior of the spacetime is defined by the gradient of the transitivity surface 
area K (t, z) = 7 a &, which can be globally timelike, spacelike, null or vary from point 
to point. We will not restrict our attention to any of these specific cases keeping the 
determinant arbitrary; just mention that in cosmology one is most interested in a globally 
timelike case, which includes unisotropic Bianchi I- VII models and their inhomogeneous 
generalizations, and a general case where the gradient of the transitivity area may vary 
from point to point. The latter includes the most general Bianchi VIII and IX types with 
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local rotational symmetry with or without inhomogeneities [12], precisely the case which 
will serve us as example in this paper. If one of the Killing vector fields £1 = d± and 
£2 = d2 is orthogonal to the hypersurface obtained by dragging the surface t-z along the 
other Killing vector field then the spacetime is globally diagonalizable. We thus start from 
such a spacetime in vacuum and notice that we can write 



lab (t, z)dx a dx b = K(t, z) e p(M) (dx 1 ) 2 + e~ p{t ' z) (dx 2 ) 2 



(2.2) 



where the function p(t, z) satisfies the linear wave equation 

f Kp 







dt dz 

and the other metric function f(t, z) may be found by a quadrature fTI ] 

One now applies one of the standard techniques to generate the non-diagonal solution 
starting from the diagonal seed ( |2.2[ ). This can be done in multiple ways: either using 
inverse scattering method of ref. | JL4|| , the solution generating algorithm described in 



or an appropriately addapted HKX transformation Ell]. Actually, even a simple Ehlers 
rotation of the Killing vectors |22| will do the job. Although this last procedure does not 
generate a genuinely non-diagonal solution (one may globally re-diagonalize the metric 
leading to a new solution), combined with a T-duality transformation in the string frame 
it will produce the desired non-vanishing components for the B^J field. 

Assuming that the non-diagonal vacuum solution has been constructed, we now have 
to solve Einstein equations for the metric (|2.1|) with a minimally coupled massless scalar 
field (our proto-dilaton) . Fortunately, these solutions may be quite easily obtained, since 
a massless scalar field has the same characteristics of propagation as gravity, and does not 
introduce any extra degree of non-linearity into the problem. In fact, if (p(t, z) is the scalar 
field, the new solution of the Einstein equations may be written as 



(2.3) 



7a6 (M) = 7a6 C (M), 

f(t,z) = r c (t,z) + r(t,z), 

where / vac (t, z) is the vacuum solution of the line element ( |2.1| ), f sc (t : z) is determined by 

2K 



r 



f 



K '2 _ K 2 

2K 
K' 2 - K 2 L 



2K 



K(<P Z + 



K (<P 2 + 



2K(f)(f) 



(2.4) 



and the scalar field 0(£, z) verifies the following linear differential equation 

0. 



(2.5) 



If the gradient of the transitivity surface area is globally timelike one may choose K ~ t, 
and the solutions of the equation may be presented as combinations of the Bessel functions 
of the first and second kind UI71 



(3 log* + C{A U cos[uo(z + z Q )]J (ujt)} + CiB^ cos[u{z + z )]N (ut)} 

' Z + Zi 



dj arc cosh 



t 



where £ indicates linear combinations of the terms in curly brackets, and u can have 
a discrete or continous spectrum. On the other hand, in the more general case when 
the gradient varies from point to point and the spatial sections have S 3 topology, as it 
happens in the case of Bianchi IX models, K ~ sin t sin z and the general solution of 
the equation (|2.5|) can be expanded in Legendre polinomials of the first and second kind 



ailog 



tan 



+ « 2 log 



tan 



+ as log | sint sin z\ 



+ J2[A e P e (cos t) +B e Q £ (cost)] [C e Pi(cosz)+DiQ e (cosz)] 



(2.6) 



e=o 



Equations (|2.3|) , ( |2.4j ) and ( |2.6|) summarize the preliminary construction in the Ein- 
stein frame. Now we transform the solution into the string frame by the conformal trans- 
formation PHjPB| 



ds* 



D 2<t>{t,z) ds < 



This provides us with a solution to four-dimensional dilaton gravity. We are, nevertheless, 
interested in having non-trivial values for other background fields, mainly the two-form 
potential B^J appearing in ( |1.1|) . To accomplish this we can use a T-duality transforma- 
tion. Taking adapted coordinates in which x° denotes the coordinate along the Killing 
vector chosen to dualize, we find new values for (G M! ,, <f>, B^J) given by |27[] 



Coo — 


1 

Goo ' 


Gon — 


Goo 


5(1) 
B 0ti ~ 


Gov 
Goo ' 


bW = 





Guv — G 



r* r* j_ rW r(i) 



>0 ^Qv 



fj,v 



G 



00 



Go^B^]} + GqvB^q 



(2.7) 



G 



oo 



log \/Go~o- 
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The new background fields so obtained automatically satisfy the field equations derived 
from the effective action ( |1 . 1|) ||27|| . Since we have two commuting Killing vectors we can 
dualize with respect to both of them to end up with non-zero values for several components 
of Bj^J. At this point it is important to stress that we are using Buscher's formulae (2.7) 
just as a formal procedure to generate new solutions of the low-energy field equations. 

So far, the application of our algorithm leads to solutions to the low-energy equa- 
tions for the "universal massless spectrum" of the four dimensional heterotic and type-II 
superstring. In the latter case we would be interested in getting solutions including also 
background values for the R-R fields. In the case of the type-IIB superstring this can 
be done using the invariance of the whole effective action S un i v + 5™ under the four- 
dimensional SL(2, R) transformation (|1 . 2| ) ||. In particular, starting with a solution in 



(2) 

which x = B^J = 0, we can generate non-trivial values for these R-R fields, thus complet- 
ing our generating algorithm. The important point here is to realize that, in this last step, 
we are doing more than just performing a S-duality transformation. Although SL(2, R) is 
a symmetry of the low-energy effective type-IIB supergravity, only a subgroup <SX(2, Z) is 
an actual symmetry of the full underlying string theory. Therefore in perfoming a generic 
SL(2, R) rotation we will end up with solutions which are not equivalent to the original 
one at the level of the string theory. 



3. Inhomogeneous Mixmaster string cosmologies 

To illustrate the use of the algorithm we suggest to look at the following example. It 
is well known in Cosmology that Bianchi type IX models (the famous Mixmaster universes 
||), let alone their inhomogeneous generalizations, represent a very interesting class of 
models to study. Due to the presence of spatial curvature and the nontrivial S 3 topology 
one may investigate their effects on the initial expansion. We consider here the locally 
rotationally symmetric (LRS) case which may be incorporated within the Einstein-Rosen 
spacetimes . To shorten the procedure we pass over the first two steps and start directly 



with the solution to Einstein-Klein-Gordon equations obtained previously by one of us JT^ . 
We can write it in the string frame as 

2M 

ds 2 = (tan0 " e 2 ^ 1 ^ ^^(-dt 2 + d9 2 ) + [I 2 (t) sin 2 9 + I 2 (t) cos 2 9} dip 2 
+ li(t)dip 2 + 2 II {t) cos 9 dip dtp} 



where the functions f(t, 9), I\{t) and I$(t) are defined by 

/(*, 9) = 21og/i(t) + C 2 sin 2 2t(l + 3 cos 2 9) sin 2 9 

+ 2C 2 sin 2 t(l - 3 cos 2 t) sin 4 0-16 —C cos t sin 2 0, 
2 ,. k 2 (tan^) ¥ + (cot|)^ 



2A (tan | + cot |) 2 
2A 



(tan |) fe +(cot|) 



2A ) 
— 



and the coordinates < < 7r, < (p < 2ir and < ip < 4ir are taken to be Euler angles. 
For the sake of simplicity, we have chosen a scalar dilaton field <p(t, 9) = (f>o(t) + (f>i(t, 9) 
containing just two modes, the homogeneous one 

, , M t 
(Po{t) = -^-logtan-, 

and the inhomogeneous growing mode (pi(t,9) compatible with the boundary conditions 
at 9 = 0,tt P], 

C 

0i(t,0) = — (l-3cos 2 t)(l-3cos 2 9). 

The S 3 topology plays an important role near the initial singularity by imposing strong re- 
strictions on the allowed modes of the scalar field, excluding for example strictly decreasing 
inhomogeous solutions. 

After T-dualizing with respect to the Killing vector £ = we trade the off-diagonal 
term in the metric, G^, for a non- vanishing component B^l of the rank-two tensor 
potential. The dual metric has finally the form 

I 2 (t) sin 2 9 + I 2 (t) cos 2 9 



d~s 2 = (tan^V * e 2 ^ \e^\-dt 2 + d9 2 ) + 



t\~ 2 ^ e -2Mt,9) 



(3.2) 



+ l tan 2 ) I 2 (t) sin 2 9 + li(t) cos 2 9 d(f 
while the new dilaton field (p(t, 9) turns out to be 



2 



-- log[/ 2 (t) sin 2 9 + li(t) cos 2 9). (3.3) 



The newly generated two-form field can be expressed as 

~ , „. I 2 (t) cos 9 

Bit, 9) = — - f-^ — — — dip A dip. 

K 1 I 2 {t) sin 2 9 + /f(t) cos 2 9 Y 



It is remarkable that the expression (|3.3| ) for the transformed dilaton is independent of the 
particular dilaton we began with. However, the latter leaves its imprint on the geometry, 
as can be seen from ( |3.2|) . 

Several interesting things happen with this solution. When the homogeneous part 
of the dilaton 0o is absent the solution is everywhere regular in the sense of curvature 
invariants and fields, whatever the value of k, except for the points located along the axis 
# = 0,7ratt = 0. This can be traced back to the fact that those points are on singular 
orbits of the Killing vector field d^. When 4>o(t) / we get a Big-Bang singularity as t 
goes to zero, as well as a Big-Crunch in t = re. With respect to the homogeneous character 
of the initial singularity, we find that for the metric before T-duality ( |3.1| ) the singularity 
is approached homogeneously, whereas (|3.2| ) gives rise to an inhomogeneous Big-Bang. As 
will be discussed later this is not unexpected after all, since some types of inhomogeneities 
are generated by the T-duality transformation. 

Suppose now the inhomogeneous part of the dilaton field is switched off, C ~ 
and (po(t) 7^ 0. The solution before T-duality describes then a Bianchi IX LRS cosmology 
which becomes inhomogeneous after the transformation. Had we started with the isotropic 
Friedman- Robertson- Walker solution (2/ 2 = 2/| = fcsint) before T-dualizing, we would 
have finished with anisotropic Kantowski-Sachs type cosmology, 

ds 2 = tan - I 2 (t)(-dt 2 + d9 2 + sin 2 6dip 2 ) 



2, 

1 / A TV ^ 

which falls into the class of solutions recently studied by Barrow and Dabrowski in |1|]. 

One thing to be learned from the above examples is that inhomogenous and/or 
anisotropic spacetimes can be related to homogeneous and isotropic backgrounds via T- 
duality. Since this is a symmetry of string theory, and therefore strings cannot distinguish 
between dual spacetimes, it seems that isotropy and homogeneity, and even some types 
of curvature singularities, become less intuitive concepts in string theory than they are 
in classical General Relativity. This is not really a surprise, since it is well known that 
string theory forces us to revise the physical role of other classical concepts, such as the 
topology of the spacetime or the value of the cosmological constant [28|. Imagine now 
an inhomogeneous spacetime T-dual to a homogeneous one, so the physics at the string 
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scale is the same in both backeroun dsQ. This seems to pose an apparent puzzle, for in the 
field theory limit (a — > 0) we are left with point particles which could, in principle, probe 
these inhomogeneities. Looking more carefully at the problem it turns out that this is 
not really the case, because of the hidden presence of the Planck length Va in Buscher's 
formulae. This makes the typical scale of these T-duality-generated inhomogeneities to be 
of the order of the Planck length, so they cannot be detected using low-energy (particle) 
probes. Thus the dual spacetime looks as homogeneous as the original one when looked 
at large distances. The moral of the story is then that there are certain types of inhomo- 
geneities, those that can be removed by a T-duality transformation, which in a sense are 
"pure gauge" from the string theory point of view. In the same vein, this generalizes to 
the case when we start with an inhomogeneous background: now we have, in addition to 
the long- wavelength inhomogeneities inherited from the original geometry, those generated 
by T-duality. Since their typical scale is of the order of the Planck length they look like 
ripples on a smooth manifold. 



4. Conclusions 

Before closing let us summarize the main results obtained. We have designed a gen- 
eral algorithm which allows the construction of cosmological solutions to the low-energy 
effective theory of the heterotic and type-IIB superstring, starting with vacuum solutions 
to Einstein equations. This provides us with a powerful tool to construct new solutions in 
string cosmology, given the huge amount of vacuum spacetimes with two or more commut- 
ing spacelike Killing vectors known in General Relativity. Starting from either homoge- 
neous or inhomogeneous seed metrics, we pass to Einstein-Klein-Gordon solutions, to end 
up after a conformal rescaling with solutions to dilaton gravity. Using T-duality and the 
S-duality of the type-IIB superstring we finally get cosmological spacetimes for both the 
heterotic and type-IIB theories, in the latter case also including non-trivial R-R fields. 

As an illustrative example we have constructed the string theory version of the in- 
homogeneous Mixmaster-type cosmology obtained in ref. We have also discussed 



the role played by inhomogeneities in string cosmology and argued that there is a class of 



7 It is important to stress that physics in dual backgrounds looks different only when described 
using momentum modes in both of them. If one, however, uses the winding states associated with 
the compact isometry in the dual background, the physical equivalence manifest itself explicitly. 
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those, the ones generated by T-duality, that have no physical relevance at any energy scale. 
Thus, T-duality classifies inhomogeneous spacetimes into physically equivalent pairs. As 
it happens with many other issues, string theory seems to challenge our physical intuition 
about the real meaning of a homogeneous geometry. 
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